Impact of r-modes on the cooling of neutron stars 
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Abstract. Studying the frequency and temperature evolution of a compact star can give us valuable 
information about the microscopic properties of the matter inside the star. In this paper we study 
the effect of dissipative reheating of a neutron star due to r-mode oscillations on its temperature 
evolution. We find that there is still an impact of an r-mode phase on the temperature long after the 
star has left the instability region and the r-mode is damped completely. With accurate temperature 
measurements it may be possible to detect this trace of a previous r-mode phase in observed pulsars. 
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INTRODUCTION 

Neutron stars are the remnants of the core collapse supernova explosions of massive 
stars which are born with temperatures higher than 10 11 K. For core temperatures larger 
than 10 6 -10 7 K the cooling of the star is dominated by neutrino emission from its interior. 
Oscillations of the star, such as r-mode oscillations, also affect the temperature evolution 
of the star by viscous heating which can drastically alter the thermal evolution. 

The r-modes are non-radial pulsations of neutron stars that are primarily driven by 
Coriolis forces and are coupled to gravitational radiation [1, 2, 3]. R-modes are damped 
by bulk and shear viscosities [4], and therefore they connect the microscopic properties 
of the matter inside the star, which depend on the low energy degrees of freedom and the 
equation of state, to the macroscopic and observable properties of the star. R-modes are 
generally defined to be the solutions of the perturbed fluid equations having (Eulerian) 
velocity perturbations of the form 

8v = aRa( K -) rxVY mm e iat + 0(Q. 3 ) (1) 

where a is a dimensionless amplitude, R and Q. are the radius and angular velocity of 
the unperturbed star, Y mm is the spherical harmonic and CO is the frequency of the mode, 
given by (O = — (m — l)(m + 2)/(m + + 0(Q. ) in the inertial frame. The r-modes 
evolve with time dependence of e^ mt ^^\ where 1/t can be decomposed as 

1 + ' + ' + ... (2) 
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Here Tg, Tb and T5 are gravitational radiation, bulk viscosity and shear viscosity 
timescales, respectively, and the dots denote other dissipative mechanisms. R-modes are 
unstable when x < 0. 



A phenomenological model describing the r-mode evolution has been given by Owen 
et. al. [5] (see also [6]). Since the amplitude of the r-mode grows exponentially in the 
unstable region, where the gravitational radiation timescale is smaller than the damping 
timescale due to viscosity, there must be some non-linear mechanisms that saturates the 
amplitude at a finite value. Supra-thermal bulk viscosity [7] is one of these non-linear 
mechanisms, but in the case that only the core of the star is considered and the effect 
of the crust is ignored, it can only saturate the r-mode at large amplitudes [8]. Mode 
coupling is another mechanism that can saturate the r-mode at small amplitude when 
hyperons and therefore large hyperon bulk viscosity is included in the computations 
[9, 10]. Core-crust boundary layer effects [11, 12] are another possible mechanism. At 
this point it is not entirely clear which mechanism is actually responsible for saturating 
the r-mode amplitude. In this paper we study the evolution of a neutron star in the 
presence of r-mode oscillations with focus on the effects on its thermal evolution. Our 
approach in studying the evolution of the star is similar to [13], but the main difference 
that distinguishes our results from previous analyses is a difference in the reheating term 
which appears in the temperature evolution. 



R-MODE EVOLUTION 

Studying the spin-down and temperature evolution of a neutron star that involves r- 
modes requires to solve a system of three coupled evolution equations for the rotation 
frequency Q. of the star and the r-mode amplitude a, as well as the temperature T. These 
equations are derived from energy and angular momentum conservation [6] 

= -2Q = - (i_a 2 <2), (3) 
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where Q = 37/ (2/) and the damping time Ty is given by 1 fxy = 1 fxs + 1/is H . The 

dots denote other possible dissipative mechanisms, like boundary layer effects. The time 
scales 1 / Tj are given by 

I = -A. (4) 

where Pq is the power radiated by gravitational waves and Pb and Ps are dissipated 
powers due to (subthermal) bulk and shear viscosity, and E m is the energy of the mode 
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TABLE 1. Radius, Kepler frequency and radial integral parameters that appear in the 
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TABLE 2. Radial integral parameters that appear in the moment of inertia, angular momen- 
tum of the mode, specific heat and neutrino emissivity of the star. 
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Here we concentrate on the lowest m = 2 r-mode whereas 5 and o are given in Table 1. 
The parameters V, S, I and J, which are given in Tables 1 and 2, encode the radial 
integration over the star and Aqcd and Aew are characteristic strong and electroweak 
scales introduced to make these quantities dimensionless. In our calculations we have 
used Aqcd = 1 GeV and Aew = 100 GeV. 

The third equation for the temperature evolution is obtained by noting that the tem- 
perature of the star decreases mainly due to neutrino emission from the interior which 
in an average mass hadronic star is dominated by modified Urea process (in a massive 
star direct Urea process are allowed in an inner core and should be included in the neu- 
trino emissivity as well), and it increases due to the dissipation of the r-mode energy by 
viscosity and other mechanisms 

where the dissipated power is P v = P$ + Pb H . Here the total (subthermal) neutrino 

luminosity L v and the integrated specific heat Cy of the star are given by 

Lv = T^^T\ C V = 4nA 3 QC lR 3 C v T v . (10) 

A EW 

The dimensionless parameters L v and Cy that involve radial integration over the star and 
and v are given in Table 2. 

Since we don't know which mechanism is actually responsible for the saturation of the 
r-mode amplitude we employ the simple model used in [5]. When solving the evolution 
equations we assume that there is a non-linear dissipative mechanism that saturates the 
mode at a fixed value a sat . The evolution is then performed in three steps. First the 
evolution equations (3) and (9) are solved until the amplitude reaches its saturation 
value a sat . Saturation requires then that there is a strong amplitude-dependent dissipative 
mechanism so that the the right hand side of the amplitude equation in eq. (3) vanishes 
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FIGURE 1. Left panel: The spindown evolution of a young 1.4 M Q APR neutron star in the angular 
frequency versus temperature plane. Different curves correspond to different saturation amplitudes, from 
10~ 4 on the left to 1 on the right. The dotted curve shows the boundary of the instability region. The 
initial conditions that we used for solving the evolution equations are Q.q = 0.8 Q.k, where £1% is the 
Kepler frequency, Oq = 10~ 6 and Tq = 30MeV. Right panel: The evolution of the angular frequency of 
the star as a function of time for the same model shown on the left panel. 



at saturation, so that 
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and the amplitude evolution stops. Therefore at saturation, in the residual evolution 
equations 1/zv nas to be replaced accordingly and the required dissipation to saturate 
the mode is entirely determined by gravitational physics. Since the dissipation that stops 
the growth of the mode inevitably heats the star, in contrast to [5] we also consistently 
replace Py by Pg/(1 ~ (X 2 Q) m me thermal equation. When viscosity alone can not 
stop the growth of the r-mode and an additional dissipative saturation mechanism is 
required this leads to a considerably stronger reheating than when only the effect of 
viscous dissipation is considered, as has been done in previous analyses [13, 15, 16]. 
The evolution in this saturation stage is performed until the boundary of the instability 
region is reached. From that point on, in a third step, the initial system of equations is 
solved to describe the subsequent decay of the r-mode and the cooling outside of the 
instability region. 

Here we study the evolution of a neutron star that is made of "APR hadron matter" 
[17] (see [14] for more details on our model). We also assume that the core of the 
star dominates the relevant quantities and neglect the effect of the crust. For our star 
model shear viscosity arises from leptonic scattering and the bulk viscosity and neutrino 
emission is dominated by modified Urea reactions. 
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FIGURE 2. Left panel: Evolution of the r-mode amplitude with time for the same model as Fig.l; right 
panel: Evolution of the temperature of the star with time. 



4.0 
3.5 

43 

1 3.0 

i 

o 1 

fcl 2.5 

QJ 

C 

J 2.0 
1.5 
1.0 





> ^^a sat =io- 4 ; 






ffsat=l 











0.0 0.5 1.0 1.5 2.0 2.5 3.0 

t/t r-mode 

FIGURE 3. Ratio of the star's temperature with and without the reheating due to r-mode dissipation 
versus time in units of the r-mode spin-down time (which increases strongly with decreasing saturation 
amplitude as seen in fig. 2). 



RESULTS AND DISCUSSION 

In Figs. 1, 2 and 3 we show the results for the evolution of a young 1.4M APR neutron 
star with the following initial conditions: = 0.8^, where Q.k = 4/9y / 27fGpo is the 
Kepler frequency, (Xq = 10 and Tq = 30MeV. We perform our calculations for differ- 
ent values of ct sat , from 10~ 4 to 1. The left panel of Fig 1 shows the rotation frequency 
of the star versus temperature. The dotted curve in that plot shows the boundary of the 
r-mode instability region for our star model where in the region above the curve the r- 



mode is unstable (Tq < Ty) and its amplitude grows exponentially until it gets saturated 
by the nonlinear mechanism, whereas in the region below the curve viscosity can damp 
the r-mode (Tg > Ty). As can be seen in Fig. 1, in the beginning of the evolution, cool- 
ing is very fast and the star enters the instability region in a fraction of a second. The 
amplitude of the r-mode grows exponentially until the amplitude reaches its saturation 
value a sat and the star reaches a thermal steady state where viscous heating equals neu- 
trino cooling. The star then spins down until it reaches the boundary of the instability 
region, where the amplitude decays. Different curves on these plots show the results of 
the evolution for different values of the saturation amplitude, from 10 to 1. Due to the 
stronger reheating the curves for larger values of a sat lie at higher temperatures. 

On the right panel of Fig. 2 the evolution of the temperature of the star versus time is 
shown. Initially the r-mode is absent and the star follows a straight line in the logarithmic 
t-T plot corresponding to mere cooling due to neutrino emission. As can be seen, at large 
saturation amplitudes a sat > 0.1 the star temporarily cools below the corresponding 
steady state curve before the amplitude is large enough so that reheating becomes 
relevant. Yet, once the r-mode is saturated the dissipative heating slows the cooling of the 
star significantly at all saturation amplitudes. After the star leaves the instability region 
and the amplitude goes to zero the star cools further by neutrino emission and approaches 
in the logarithmic plot the initial linear behavior, so that the r-mode effectively merely 
delays the cooling evolution. The effect of reheating is even larger for smaller values of 
Ct sat because in those cases the spin-down takes longer and the star spends more time in 
the unstable region, increasing the delay effect. 

The dissipation due to bulk viscosity eq. (6) and the neutrino luminosity eq. (10) were 
given in the subthermal approximation where non-linear amplitude-dependent terms are 
neglected. Taking into account the suprathermal non-linear amplitude-dependence [7] 
hardly affects the dissipation when there is another non-linear mechanism that saturates 
the mode due to even stronger dissipation. Although less obvious, we also find that 
the amplitude dependence of the neutrino luminosity [18, 19] hardly affects the results 
since the strong reheating due to the dissipative saturation mechanism strongly delays 
the cooling and keeps the star at larger temperatures where suprathermal effects are 
small. This point will be discussed in more detail elsewhere. 

To analyze the reheating effect due to the r-modes in detail we compare in Fig. 3 the 
thermal evolution with and without r-mode reheating. Shown is the ratio r r _ mo d e over 
r no _ r _ mode , where r nor _ mode describes the cooling in the absence of r-mode oscillations, 
versus time in units of the r-mode spin-down time ? r - mode* which is the time when 
the star leaves the instability region. The plot shows that the reheating increases the 
temperature during the r-mode phase by a constant factor which can take values up to 
5 for realistic amplitudes. Whereas for large amplitudes the evolution in the vicinity 
of the instability region is affected by the simplified model prescription to keep the 
amplitude constant until the boundary is reached (where it should actually vanish), for 
sufficiently low amplitudes 0t sat < 0.1 the thermal evolution outside of the instability 
region takes in the scaled coordinates of Fig. 3 a universal form that is independent of 
the amplitude. The temperature decays with time scales that are considerably longer 
than the decay of the r-mode amplitude so that effects of the previous r-mode phase 
on the thermal evolution are still present long after the r-mode oscillation ceased. This 
is interesting since, nearly all observed young pulsars spin already too slow that they 



could be unstable to r-mode oscillations at present, as discussed in more detail in a 
forthcoming publication [20]. There we show that r-modes are a consistent mechanism 
to explain these low frequencies for saturation amplitudes 0.001 < 0C sat < 0.1 for which 
r-modes are present over a significant fraction of the age of young stars. Fig. 3 shows 
that even after twice the r-mode spin-down time, the temperature of the star that had 
an r-mode history is still more than ten percent larger than that of a star without. This 
presents a unique signature of a previous r-mode phase that is nearly independent of the 
unknown saturation amplitude. Although this effect is at present overshadowed by larger 
observational and model uncertainties, more robust data and a better understanding of 
these systems could in the future allow to detect a past r-mode phase in temperature 
measurements of young pulsars. 
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